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Abstract 

We describe the physical hypotheses underlying the derivation of an approximate model of 
water waves. For unidirectional surface shallow water waves moving over an irrotational 
flow as well as over a non-zero vorticity flow, we derive the Camassa-Holm equation by an 
interplay of variational methods and small-parameter expansions. 
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1 Introduction 

The Camassa-Holm equation reads 

Ut + 2kU^ + 3UU^ - Ut^^ = 2U^U^^ + UU^^^ (1) 

with x e R, t e R, [/(a;, i) G R. Subscripts here, and later, denote partial deriva- 
tives. The constant k is related to the critical shallow water speed. For k—0 the 
equation ([T]) possesses peaked soliton solutions [S]. The physical derivation of ([T|) 
as a model for the evolution of a shallow water layer under the influence of grav- 
ity, is due to Camassa and Holm [5]. See also Refs. [5], [TH], [53] for alternative 
derivations within the shallow water regime. In [31j it was shown that the equation 
([T]) describes a geodesic flow on the one dimensional central extension of the group 
of diffeomorphisms of the circle (for the case k=0, a geodesic flow on the diffeo- 
morphism group of the circle, see also [H]). It should be mentioned that, prior to 
Camassa and Holm, Fokas and Fuchssteiner [SOj obtained formally, by the method 
of recursion operators in the context of hereditary symmetries, families of integrable 
equations containing ([T]) . These equations are bi-Hamiltonian generalizations of the 
KdV equation and possess infinitely many conserved quantities in involution |20| . 
In order to see how the equation ([T]) relates with one of equations in the families 
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introduced in [50], see [52], §2.2. Also, the equation (5.3) in [3T] is the equation ([T|) 
but with errors in the coefficients. 

The Camassa-Holm equation attracted a lot of interest, due to its complete 
integrability 0, [15], [25, [33] (for the periodic case), [2], [2], [E] and the citations 
therein (for the integrability on the line), the existence of waves that exist for all 
times as well as of breaking waves |1JJ and the presence of peakon solutions [5]. 
Of interest is also the fact that its solitary waves are solitons, recovering their 
form and speed after nonlinear interaction [3], |28| . These waves are smooth if 
K 7^ and peaked if k = 0, with both wave forms being stable and thus physically 
recognizable |16j , [17] . An important question to be answered is whether solutions 
of the water wave problem can really be approximated by solutions of the Camassa- 
Holm equation. In ^0 it is shown that suitable solutions of the water wave problem 
and solutions of the Camassa-Holm equation stay close together for long times, in 
the case e — S^, e being the amplitude parameter and S the shallowness parameter. 
In [30j it is also shown that the peakon equation cannot strictly be derived from 
the Euler equation and hence it is at most a phenomcnological model. 

In the references [TH] , [H] , [H] , the derivation of the Camassa-Holm equation is 
made under the assumption that the fluid motion is irrotational. This assumption 
does not appear explicitly in the initial derivation [5]. This derivation is based on 
the Green-Naghdi equations and on the choice of a submanifold in the Hamiltonian 
representation of the Green-Naghdi equations. For the Green-Naghdi equations, one 
does not have to impose any condition that the fluid is irrotational, but, as shown 
in [26] , in order to have a consistency for the second assumption one should impose 
an irrotational condition on the flow. Thus, another important question to be 
answered is the question whether a Camassa-Holm equation can be obtained when 
the fluid below the free surface has non-zero vorticity. By examining the underlying 
flow one has some important and new consequences. In [27] it is considered the 
Camassa-Holm equation for water waves moving over a shear flow. The calculation 
is completed for the case of a linear shear, that is, the underlying flow has constant 
vorticity. In this case it is shown that the equation for the surface wave is not a 
Camassa-Holm equation, but the Camassa-Holm equation can exist for a simple 
nonlinear function of the horizontal velocity component of the perturbed flow field, 
at a certain depth. 

The purpose of the present paper is to derive by a variational approach in 
the Lagrangian formalism, the Camassa-Holm equation for unidirectional surface 
shallow water waves when the fluid below the free surface has non-zero vorticity. 
We will consider in turn the cases of an irrotational flow (for this case see also 
[53]), a rotational flow with constant vorticity and finally, an arbitrary flow. For 
the linear shear case, we show that the displacement of the free surface from the 
flat state, satisfies the Camassa-Holm equation ([Ij with k ^ 0. For an arbitrary 
unidirectional shallow water flow, we get that the displacement of the free surface 
from the flat state, satisfies a generalized Camassa-Holm equation. 

2 Nondimensionalisation and scaling of the governing equations 
for water waves 

We consider water moving in a domain with a free upper surface at z — ho + ri{x, t), 
for a constant ho > 0, and a flat bottom at z = 0. The undisturbed water surface is 
z = ho. Let {u{x , z , t) , v{x , z , t)) be the velocity of the water - no motion takes place 
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in the y-direction. The fluid is acted on only by the acceleration of gravity g, and 
the effects of surface tension are ignored. For gravity water waves, the appropriate 
equations of motion are Euler's equations (EE )(see [25]). Another realistic as- 
sumption for gravity water wave problem is the incompressibility (constant density 
p) (see [29]), which implies the equation of mass conservation (MC). The bound- 
ary conditions for the water wave problem are the kinematic boundary conditions 
as well as the dynamic boundary condition. The kinematic boundary conditions 
(KBC) express the fact that the same particles always form the free water surface 
and that the fluid is assumed to be bounded below by a hard horizontal bed z — Q. 
The dynamic boundary condition (DBC) expresses the fact that on the free surface 
the pressure is equal to the constant atmospheric pressure denoted Pq. Summing 
up, the exact solution for the water-wave problem is given by the system 

Ut + uu^ + vu^ = -^Px 
Vt + UVx + VVz = -^Pz - 9 

Ux + Vz = 

V = Tjt + urjx on z — + ri{x, t) 
V = Q on z = 



(MC) 
(KBC) 



(2) 



p = po, on z — ho + r^{x,t) (DBC) 

where p{x^ z, t) denotes the pressure. While a mathematical study of exact solutions 
to the governing equations ^ for water waves can be pursued (see for example [32] 
for the periodic steady solutions in the irrotational case and |18| for the periodic 
steady solutions in the case of non-zero vorticity), to reach detailed information 
about qualitative features of water waves it is useful to derive approximate models 
which are more amenable to an in-depth analysis. 

We non-dimensionalise the set of equations ^ using the undisturbed depth 
of water h^, as the vertical scale, a typical wavelength A, as the horizontal scale, 
and a typical amplitude of the surface wave a (for more details see [25], [M])- An 
appropriate choice for the scale of the horizontal component of the velocity is ^/gho- 
Then, the corresponding time scale is and the scale for the vertical component 

of the velocity is ho^^. Thus, we define the set of non-dimensional variables 




/gho 



gho"' (3) 



where, to avoid new notations, we have used the same symbols for the non-dimensional 
variables x, z, rj, t, u, v, in the right-hand side. The partial derivatives will be re- 
placed by 



qhn \/ Qh() \/ Qhn 

Ut ^ ^Ut, Ux ^ ^^Ux, Uz l-> ^%f-Uz, 



Let us now define the non-dimensional pressure. If the water would be stationary, 
that is, u = V = 0, from the first two equations and the last condition with 77 = 0, 
of the system ©, we get for a non-dimensionalised z, the hydrostatic pressure 
Po + P5^o(l ~ z). Thus, the non-dimensional pressure is defined by 

P Po + pghoil - z) + pghop (5) 
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and 

Px ^ P~P^' '"^ + P9P^ (6) 

Taking into account ([3]) , ([4]) , ([5]) , and ([6]) , the water-wave problem ([2]) writes in 
non-dimensional variables, as 

Ut + UUx + VUz = -px 
S'^ivt + UVx +VVz) ^ -Pz 

Ux + Vz = (7) 
V = e{r]t + urjx) and p = erj on z = 1 -|- er]{x, t) 

V = Q on z = 

where we have introduced the amplitude parameter ^ — and the shallowness 
parameter 6 — The small-amplitude shallow water is obtained in the limits 
e — » 0, 5 — > 0. We observe that, on z = 1 -t- ery, both v and p are proportional to e. 
This is consistent with the fact that as e — ^ we must have v ^ Q and p ^ 0, and 
it leads to the following scaling of the non-dimensional variables 

p I— > ep, {u, v) 1-^ e{u, v) (8) 

where we avoided again the introduction of a new notation. The problem ([7]) be- 
comes 

Ut + e{uux + vuz) = ~Px 
S^[vt + e{uvx + vvz)] = -pz 

Ux + Vz = Q (9) 
V — "qt + (-urix and p — rj onz = l + trj{x^ t) 

V — Q on z 

The two important parameters e and 5 that arise in water-waves theories, are used 
to define various approximations of the governing equations and the boundary con- 
ditions. The amplitude parameter e is associated with the nonlinearity of the wave, 
so that small e implies a nearly-linear wave theory. The shallowness parameter 5 is 
associated with the dispersion of the wave, it measures the deviation of the pressure, 
in the water below the wave, away from the hydrostatic pressure distribution. 



3 Variational derivation of the Camassa-Holm shallow water 
equation 

In contrast to the KdV equation, which is a classical integrable model for shallow 
water waves, the Camassa-Holm equation possesses not only solutions that are 
global in time but models also wave breaking. The only way that singularities can 
arise in finite time in a smooth solution is in the form of breaking waves, that is, 
the solution remains bounded but its slope becomes unbounded [7] . Even if a wave 
breaks, there is a procedure to continue uniquely the solution after wave breaking 

SI- 

There are different approaches which lead to the Camassa-Holm equation in 
the shallow water regime. The original derivation of this equation [5| consists in 
making approximations by relating m := u — Uxx and rj in the Green-Naghdi Hamil- 
tonian system and preserving the momentum part of its Hamiltonian structure. The 
Camassa-Holm equation can be also obtained from the governing equations ([9]) , by 
the use of a double asymptotic expansion, valid for e — > 0, (5 — > 0, retaining terms 
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0(1), 0(e), 0((5^), 0{e6^) (see for example [55], [H]). As a result, a single nonlinear 
equation for rj will be obtained and thus all variables will be expressed through the 
solution of this equation. In what follows we will consider a derivation of the model 
from the governing equations for the water-wave problem by a variational approach 
in the Lagrangian formalism [8] , [2l| ■ 

One observes that the parameter 6 can be removed from the system © (see 
pS]). this being equivalent to use only ho as the length scale of the problem. In 
order to do this, the non-dimensional variables x, t and v from ([3]) are replaced by 

a/c -v/e 6 

X 1-^- — ^a;, t i-^ — i, V i— > —pV (10) 
6 y/e 

Therefore the equations in the system ([9|) are recovered, but with 5"^ replaced by e 
in the second equation of the system, that is, this equation writes as 

e[vt + e{uvx + vv^)\ = -p^ (11) 

The classical approximation is the linearized problem obtained by requiring the 
amplitude of the surface to be small letting e ^ 0. The system ([9]) with the second 
equation given by ()11|) becomes linear: 

Mt + Px = 

u, + v,^Q (12) 
V = r]t and p ~ rj on z = 1 
V — on z = 

3.1 The case of an irrotational flow 

Firstly, we consider that the fluid flow is irrotational, i.e. the flow has zero vorticity. 
Then, in addition to the system ([2]), we also have the equation 

Uz-Vj;^0 (13) 

Here the velocity components u and v are written in the physical (dimensional) 
variables. If we non-dimensionalise this equation using ([3]), ([4]), we obtain 

Uz = S^v^ (14) 

After scaling ([5]) and transformation pop . the equation writes as 

Uz = evx (15) 

Therefore, in the limit e ^ 0, we get in addition to the system p^ . the equation 

Uz = (16) 

Thus, under the assumption that the fluid is irrotational, u is independent of z, 
that is, 

u — u(x, t) (17) 

From the second equation in (fT2|) . we get that p also does not depend on z. Because 
p — rj{x, t) on z = 1, we have 

p = ?7(x, t) for any < z < 1 (18) 
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Therefore, using the first equation in p2p . and taking into account (fT7|) . we obtain 

'q^{x,t)dt + T{x) (19) 



where J- is an arbitrary function. Differentiating (|19p with respect to x and using 
the third equation in (|12p . we get, after an integration against z, 



(^J Vxx{x,t)dt-P{x)^ (20) 



We underhne the fact that in our approximation the vertical velocity component 
maintains a dependence on the z-variable. In analyzing the motion of the fluid 
particles, this means that the particles below the surface may perform a vertical 
motion. This is in agreement with recent general results obtained in [10. for the 
periodic steady waves (Stokes waves) and in [12] for the solitary waves, both being 
solutions of the full Euler equations. 

Making z = 1 in (j20p . and taking into account that w = 77^ on z = 1, we get after a 
differentiation with respect to t, that rj has to satisfy the equation 

Vtt -rjxx = (21) 

The general solution of this equation is i]{x, t) — f{x — t) + g{x + t), where / and 
g are differentiable functions. It is convenient first to restrict ourselves to waves 
which propagate in only one direction, thus, we choose 

ij{x,t)^ f{x~t) (22) 

Therefore, for u and v in (fT9|) . (|20|l we have u{x,t) = r]{x,t) + J^{x), v{x,z,t) = 
—z (rjx + !F'{x)). The condition t; = yyt on z = 1, yields 

T{x) = const :— co (23) 

Thus, for the irrotational case the solution of the system (|12p plus the equation 
(|16p . can be written into the form 

T]{x,t) = f{x -t), u{x,t) = r]{x,t) + Co, v{x,z,t) ^ -zi]x{x,t) ^ ~zu^, (24) 

The solutions to the shallow water problem are determined by the evolution of 
the function rj{x,t), which represents the displacement of the free surface from 
the undisturbed (flat) state. The solution ((24|l describes the linear, non-dispersive 
surface wave. 

By consistently neglecting the e contribution, we will derive using variational 
methods in the Lagrangian formalism (see [S], [21]), the equation ^ governing 
unidirectional propagation of shallow water waves. 

In the Lagrangian picture of a mechanical system, one focuses the attention on 
the motion of each individual particle of the mechanical system. We denote by M 
the ambient space whose points are supposed to represent the particles at i = 0. 
A diffeomorphism of M represents the rearrangement of the particles with respect 
to their initial positions. The set of all diffeomorphisms, denoted Diff(M), can be 
regarded (at least formally) as a Lie group. The motion of the mechanical system 
is described by a time-dependent family of orientation-preserving diffeomorphisms 
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7(i, •) G Diff(M). For a particle initially located at x, the velocity at time t is 
7t(t,x), this being the material velocity used in the Lagrangian description. The 
spatial velocity, used in the Eulcrian description, is the flow velocity w{t, X) — 
7t(i,x) at the location X — j{t,x) at time t, that is, w{t,-) — ° 1^^- In the 
Lagrangian description, the velocity phase space is the tangent bundle TDifl[(M). 
In the Eulerian description, the spatial velocity is in the tangent space at the identity 
Id of DifF(Af), that is, it is an element of the Lie algebra of Diff(M). The Lagrangian 
£ is a scalar function defined on TDiff(il/) and the equation of motion is the equation 
satisfied by a critical point of the action a{j) — ^{j, Jt)dt defined on all paths 
{7(t, •), t G [0, T]} in Diff(M) having fixed endpoints. 

For our problem, we saw in ([M)) that the vertical component v of the velocity 
is completely determined by the horizontal component u{x, t) of the velocity, which 
is a vector field on R, that is, it belongs to the Lie algebra of Diff(R). Thus, in 
the Lagrangian formalism of our problem we take M = R and add the technical 
assumption that the smooth functions defined on R with value in R vanish rapidly 
at ±cxD together with as many derivatives as necessary (see [7j for a possible choice 
of weighted Sobolev spaces). For the one-dimensional periodic motion, one takes 
M = the unit circle. In what follows we focus on the first situation. 
Following Arnold's approach to Euler equations on diffeomorphism groups [1], the 
action for our problem will be obtained by transporting the kinetic energy to all 
tangent spaces of Diff(R) by means of right translations. For small surface eleva- 
tion, the potential energy is negligible compared to the kinetic energy. Taking into 
account (|24)) . the kinetic energy on the surface is 



-1 /"OC -1 /"OO -1 /"OO 

K=\ {u^ + v')dx = \ [u' + {l + e^ful] dx^\ {u' + ul) dx 

^ J —oo ^ J —oo ^ J — oo 

(25) 

to the order of our approximation (see [8]). We observe that if we replace the path 
7(t, ■) by 7(t, •) o i/'(-), for a fixed time-independent ij} in Diff(R), then the spatial 
velocity is unchanged 74 o 7^^. Transforming K to a right-invariant Lagrangian, 
the action on a path 7(i, •), t e [0, T], in Diff(R) is 

a(7) = 7^ / {(7* o 7-')' + mit o -/~')r}dxdt (26) 

The critical points of the action (I^S)) in the space of paths with fixed endpoints, 
verify 

-fa(7 + e^) =0, (27) 

de 6=0 

for every path (p{t, •), t g [0, T], in Diff(R) with endpoints at zero, that is, 1^9(0, •) = 
= (p{T, •) and such that 7 -I- £(/? is a small variation of 7 on Diff(R). Taking into 
account (|26p . the condition ((27|l becomes 



(7* o 7-^) I 



e=0 



[{-ft + sipt) ° h + ^] 



+5,(7*07-!)-^ [d, {{jt + eift) o + Eip)-^)] [ dxdt ^ (28) 
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After calculation (see [5]), we get 
d 



de 



e=0 



[{jt + £Vt) ° h + ^v) ^]^Vt°l ^-{ipoj ^)dx{'^t°l ^) 

-((^07-1)9,(7407-!) (29) 



and 



de 



e=0 



[a, ((7t + e(p)o (7 + e^) i)] 



dx{vt°l ^)-dx{-ftoj i)9,((po7 1) 

-(<^o 7-1)5^(7, o7~i) 

dtx{^ o 7-I) + (74 o j-^)dl{ip o 7-I) 



where are used the formulas of the type 



d 
de 



e=0 



Itx o 7 
7x o T 



Ix o 7~ 
1 



(30) 



(31) 



(32) 



9t((y5 07 i) = (^t0 7 i + (<^:c0 7 i)9f(7 i)=^i07 1- (7*0 7)^3; ((y5 07 1) (33) 
Thus, denoting 74 o 7-^ = u, from (p9)l . (jSO]) . the condition (p8)) writes as 

/ [u [dt{ip o 7"i) + udx{^ o 7"i) - ((^ o 7"l)U:r] 

J —00 

+Ux [dtxiv o 7"^) + ud'^i'P ° 7"^) 7^^)uxx] } = (34) 

We integrate by parts with respect to t and x in the above formula, we take into 
account that ip has endpoints at zero, the smooth functions defined on R with values 
in R, together with as many derivatives as necessary, vanish rapidly at ±00, and 
we obtain 



T poo 



((^07 1) [Wt + 3uUx - Utxx - '^UxUx 



c] dxdt = 



(35) 



Therefore, we proved: 



Theorem 3.1: For an irrotational unidirectional shallow water flow, the horizontal 
velocity component of the water u(x,t) satisfies the Camassa-Holm equation (Q]) for 

K = 0. 

Let us see now which equation fulfill the displacement r]{x, t) of the free surface 
from the flat state. Taking into account (IM)) . the kinetic energy on the surface is 



K 



{u^ + v^)dx 



[iv + Co)' + (1 + ev)^7i] dx 



iv + co) +vl 



dx 



(36) 
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to the order of our approximation. Transforming i^T to a right-invariant Lagrangian, 
the action on a path r{t, •), t € [0,T], in Diff(R) is 

a(r) ^ \ r r {{Tt o + co)2 + [d^{Tt o T-^)Y)dxdt (37) 

The critical points of the action (I37p in the space of paths with fixed endpoints, 
verify 

-^a(r + e$) 



de 



= 0, 



(38) 



for every path <i>(i, •), t e [0,T], in Difi'(R) with endpoints at zero, that is, $(0, •) — 
= •) and such that F + e$ is a small variation of F on Diff'(R). Taking into 
account (|37p . the condition ((38|l becomes 



fF*oF-i 



-9,(rtor- 



CO 



:=0 



[{Tt + e^t) o {T + e^y^] 

=0 

[9, ((Ft + e^t) o (F + e$)"i)] | dxdt = (39) 



After calculation, denoting Ff o F ^ — r], ([55)1 writes as 

{(^7 + Co) [dt{<i> o F-i) + 779,($ o F-i) - ($ o F-i)?7,] 

[at.($ o F-i) + 77a2($ o F-i) - ($ o F-i)7],,] } dxdt = (40) 

We integrate by parts with respect to t and x in the above formula, we take into 
account that $ has endpoints at zero, the smooth functions defined on R with 
values in R, together with as many derivatives as necessary, vanish rapidly at ±oo, 
and we obtain 

($ o F~^) [r]t + irji^x + 2co772: - rjtxx - '^VxVxx - VVxxx] dxdt = (41) 

/O J -oo 

We thus have: 

Theorem 3.2: For an irrotational unidirectional shallow water flow, the displacement 
r]{x, t) of the free surface from the flat state, satisfies the Camassa-Holm equation 
^ for K = cq. 



3.2 The case of a linear shear flow 

We assume now that the underling flow is rotational with a constant vorticity luq, 
that is, we are in the case of linear shear flow. Then, in addition to the system ([2]), 
we also have the equation 

Uz — Vx ^ const :— ujq (42) 

If we non-dimensionalise this equation using ([3]), we scale using ([5]) and we 
transform using (|10p . the equation writes as 



, '^o\/po 

Uz = eux H 

9 



(43) 
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Therefore, for the case of a Unear shear flow, in the hmit e ^ 0, we get in addition 
to the system (fT2)) . instead of (fTH)) . the equation 



9 

The relation ^5)1 remains the same but instead of (|19p we have now 



(44) 



rix[x, t)dt + J- (x) H z 

9 



(45) 



where T is an arbitrary function. Using the third equation in p^ . we get again 
([^0]) . Following the same procedure as in the irrotational case, presented after the 
relation (|20p . we obtain the solution of the system plus the equation ([44]) into 
the form 

r]{x,t) — f{x ~ t), u{x, z,t) — r]{x,t) -\ ° 2 + cq, v{x, z,t) = —zrixix,t) 

9 

(46) 

We will derive using the variational methods in Lagrangian formalism, the equa- 
tion of the water's free surface r]{x,t). Taking into account (|^ . we observe that 
ri{x,t) determined completely the velocity components u and ti. r]{x,t) can be re- 
garded vector field on R, that is, it belongs to the Lie algebra of Diff(R). 
The kinetic energy on the surface is 



JC = 



{u^ + v^)dx = 



[V + ^^^^{1 + e^) + cor + (1 + evrvl } dx 



{V + +co) +Vx 

9 



dx 



(47) 



to the order of our approximation. Transforming /C to a right-invariant Lagrangian, 
the action on a path r{t, •), t € [0,T], in Diff(R) is 



^ JO J-oo 



-1 , t^oVff^ 



+ cof + [dx{Ttor-^)r}dxdt (48) 



The critical points of the action ([48]) in the space of paths with fixed endpoints, 
verify ([38|) . for every path $(i, •), t E [0,T], in Diff(R) with endpoints at zero. 
Taking into account (|48|) . the condition ([38]) becomes 



T 

^-oo 



Ft o r H h Co 



d 

de 



e=0 



[(Ft+e$t)°(r + e$)~'] 



+5^(Ft o F-i)-^ [9, ((Ft + e$) o (F e*)-^)] 1 dxdt = (49) 
de E=o J 



After calculation, dcnotiEgo F ^ = i], writes as 



J-oo 



5 



-f Co) [c>t($ o F-1) + 77a,($ o F-i) - ($ o r-')r,x] 



+Vx [9te($ o F-i) + r/a2($ o F-i) - ($ o F-i)r/,,] } dxdt = (50) 
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We integrate by parts with respect to t and x in the above formula, we take into 
account that $ has endpoints at zero, the smooth functions defined on R with 
values in R, together with as many derivatives as necessary, vanish rapidly at ±00, 
and we obtain 



Jo J -00 



,0 ,01 ^oVPo , 1 
Vt + Srjrjx + 2 ( \-Co ] Tlx- rjtxx 



'^rjxVxx - rnxxx] dxdt = (51) 



Thus, we get: 



Theorem 3.3: For a rotational unidirectional shallow water flow with constant vor- 
ticity, the displacement rj{x, t) of the free surface from the flat state, satisfies the 
Camassa-Holm equation for k — ""^9^° -|- cg. 

3.3 The case of an arbitrary flow 

We consider now an arbitrary underlying flow. From the second equation in (|12p . 
and from p = ri(x,t) on z = 1, we get (jl8l) . Therefore, using the first equation in 
P^ . we obtain 

r]^{x,t)dt + J'{x,z) (52) 



where is an arbitrary function. Differentiating (|52p with respect to a; and using 
the third equation in (jl2p . we get, after an integration against z, 

r]xx{x,t)dt - g{x, z) + g{x,0) (53) 

where Qz{x,z) = J^x(x,z) and we have also taken into account the last condition 
in the system ([T^ . Making z = 1 in ([55)1 . and taking into account that v = rjt 
on z = 1, we get after a differentiation with respect to t, that rj has to satisfy the 
equation (HI]). We restrict ourselves to waves which propagate in only one direction, 
thus, we choose 77 into the form ((22)) . Therefore, for u and v in ([52| . ([SS)) we get 

u{x, z,t) ~ rj{x,t) + T{x, z), v{x, z,t) ^ —zrjx{x,t) — Q{x, z) + Q{x,0) (54) 

with 

g,{x,z) ^T4x,z), gix,i) = g{x,o) (55) 

arbitrary functions. 

Using the variational methods in Lagrangian formalism, we will derive now the 
equation of the water's free surface r]{x, t), which can be regarded as a a vector field 
on R, that is, it belongs to the Lie algebra of Diff(R). Taking into account ([5^ . 
and the second condition in (|55[) . the kinetic energy on the surface will be now 



T poo ^ f°° 

JC = - iu' + v')dx = - {[ii + T{x,{l + em^ 

^ J — oo ^ J -oo 

+ [-(1 + £77)77, - g[x, (1 + er;)) + g{x, Q)f]dx 

1 1'°° r 1 
^ 2l_ + + Vl}dx (56) 
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to the order of our approximation. Transforming /C to a right-invariant Lagrangian, 
the action on a path T{t, ■), t £ [0, T], in Diff(R) is 

a{T) = -j^ j |[r,or ^T{x,\)\ +[9,(r,or )fyxdt (57) 

The critical points of the action (|57p in the space of paths with fixed endpoints, 
verify ((38)) . for every path •), t G [0,T], in Diff(R) with endpoints at zero. 
Taking into account (l57|) . the condition ([38]) becomes 



T poo 



TtoT ' + j^{x,i] 



d 




de 


e=0 - 



+9.(r,or"')A 

d£ 



After calculation, denoting j. t 

i-T 



e=0 



dxdt 



d,({Vt + e^)o{V + e^)-^ 
= 77, ([5S|) writes as 

{[?/ + T{x, 1)] [dt{^ o r^') + r^d^i^ o r^') - ($ o r"')77. 

+r;, o r^') + ridli^ o r^') - o r~')77,,] } dxdt - 0(59) 

We integrate by parts with respect to t and x in the above formula, we take into 
account that $ has endpoints at zero, the smooth functions defined on R with 
values in R, together with as many derivatives as necessary, vanish rapidly at ±00, 
and we obtain 

- / (^oF )[rit + ?,riJi,^^2T{x,\)ri,^^ T.:,{x,\)r] 

Vtxx 

Jo J -00 

-'^VxTlxx - mxxx] dxdt = (60) 

In conclusion, we proved: 

Theorem 3.4: For an arbitrary unidirectional shallow water flow, the displacement 
r]{x,t) of the free surface from the flat state, satisfies a generalized Camassa-Holm 
equation 

F\x)U + Ut + [3C/ + 2F{x)] - Utxx - 2C/,C/,, + UU^^^ (61) 
with a; G R, t e R, U{x, t) e R, F{x) £ R. 
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